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PREFACE 


EXCEDE II is an Air Force Geophysics (AFGL) research program to study 
auroral infrared emission processes using a rocketbome electron accelerator 
to generate artificial aurorae. The effort is sponsored by Defense Nuclear 
Agency (DNA). The EXCEDE concept was originally proposed by Dr. A. T. Stair 
of AFGL in 1971 and since that time four successful experiments have been 


conducted. These experiments have been directed by Mr. hobert O'Neil also 


of AFGL. 
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INTRODUCTION 


EXCEDE II is the fourth in a series of research experiments con- 
ducted by Air Force Geophysics Laboratory (AFGL) to study the infra- 
red emission processes induced by rocketborne electron accelerators. 
The excitation of the upper atmosphere by the electron beam simulates 
aurora induced processes but provides the advantage that the dosing is 
more precisely quantified than in the natural auroral case. In addition, 
optical/infrared instruments carried on-board the same payload with 
the accelerator provide controlled viewing geometry of the induced 
emissions and measurements do not suffer from long path atmospheric 
attenuation as would be the case for some ground based observations. 
The other EXCEDE experiment programs were PRECEDE, EXCEDE II Test, 
and EXCEDE: SWIR (O'Neil 1976). 


The EXCEDE II payload incorporated four electron guns programmed to 
function in different combinations and capable of producing a 10 ampere 
electron beam when activated simultaneously. The optical/infrared 
instruments on-board the payload were primarily for measuring spectral data 
and included; liquid nitrogen (LN,) interferometer, LN, circular 
variable filter (CVF) spectrometer, liquid helium (LHe) CVF spectrometer, 
IN, radiometer, ultraviolet/visible spectrometer, a camera and five 
photometer units. In addition the payload contained several ionospheric 
probes and diagnostic sensors to measure ion and electron densities and 
fields. The optical/infrared instruments viewed the beam region through- 


out the flight at a stablized zenith angle to local vertical of 45°. 


The EXCEDE II payload was launched from Pad 4, Poker Flat Research 
Range, Alaska on 25 October 1979 at 0546:40 UT hours. The payload reached a 
~128 km apogee at ~187 seconds after launch. Although some difficulties 
were encountered with electron gun arcing the experiment was successful and 
data were obtained for gun on and off conditions throughout the flight 
(O'Neil 1980). 


The primary purpose of this report is to describe in detail the 


AE ee EEG ng SD 


techniques used to process and analyze the data measured with the two CVF 
spectrometers on the EXCEDE II payload. Part A of the report deals with 
the LN, CVF spectrometer (designated NS-6-4), Part B describes techniques 
employed with the LH CVF spectrometer (designated HS-3B-1) and Part C 
presents recommendations for final data processing of the HS-3B-1 CVF. 


Processed data for this instrument presented in this report is preliminary. 


? 
- 


2.0 


PART A 


IN, CVF SPECTROMETER 


Description of Nitrogen CVF(NS-6-4). 


Construction details and method of operation of the rocketborne 
CVF spectrometer have been reported by Rogers, 1976, and Wyatt, 1978. 
The CVF instruments were fabricated by the Electrodynamics Laboratory 
of Utah State University. In brief the monochromator portion of the 
sensor consists of two 180-degree circular variable (wedge) filters 
mounted so that the bandpass interference layers rotate just above the 
detector field stop which performs the function of a spectrometer 
slit (see Figure 2.0-1). The spectrometer scans a continuous wave- 
length region depending on which half of the filter is rotating over the 
detector. Spectral resolution (AA) depends on slit width and f number. 
Pertinent figures of merit for the EXCEDE NS-6-4 instrument are given 
in Table 2.0-1. 


Two important parameters for CVF data processing are the start 
times of the beginning of each spectral scan and the speed of rotation 
of the CVF. Two assumptions made are that the speed of rotation is 
constant throughout the flight and that the wavelength of radiation 
being measured is linearly proportional to time (or filter position). 
One consequence of this CVF mode of operation is that since the payload 
is constantly changing altitude the spectra obtained are associated with 
different altitudes. Each part of a spectral scan is for a different 
altitude and thus altitude must be computed for every spectral data 


point. 
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Table 2.0-1 


PERFORMANCE PARAMETERS 
FOR CVF NS-6-4 


ee aan a le Rl caer a a 


Wavelength Coverage (um) 2.1 to 3.5 , 3.25 to 35.4 
Spectral Resolution (ym) 0.05 € 2.7 um, 0.08 ¢ 4.3 um 
Detector Type InSb 

NESR (w-em™“-sr7!-um-4) 2.4 « 19710” 

Dynamic Range 2» 10° . 

Half Angle Field of View (deg) 1.871 

Spectral Scan Rate ~2/sec 


* 
reference 5 


sk 


NS-6-4 Calibration Parameters 


Radiance and wavelength calibrations for the NS-6-4 CVF were 
accomplished by both AFGL and Utah State University (Wyath 1978) prior to 
rocket launch. These calibrations showed good agreement so those of AFGL 


were used for this work. 


2.1.1 Wavelength Calibration 


The NS-6-4 CVF had two halves designated ‘ong (subscript 1) 
and short (designated s) which refers to long wavelength and short 
wavelength. The AFGL calibration equations were adjusted using the 
data and known emission features such as the CO, (v3) peak at 4.28um. 


The appropriate equations are 


(8) = 1.693 + 6.916 x 107) Py 


(2-1) 
6.79 CP. £ 46.37 , 
Ee = Bx 100 (2-2) 


Where n = number of the data point in the scan 
N = Total number of data points in the scan 


P,= percentage of n data point of total scan 
For the long half of the filter. 
a 
Ane)" -0.292 + 6.860 x 10 in (2-3) 
59.42 SP 4 97.82 
In the NS-6-4 the short half of the filter was scanned first. 


2.1.2 Radiance Calibration 


Radiance computations were done in units of megarayleighs per 
micron (MR/um) or brightness (B) since these units were used in 
calibration and permit direct physical interpretation of the results, 
In some cases the term "'existence' is also used (Wyath 1978). To 
convert telemetry volts Vay contained in the data channels the 


following calibration equations apply 


" . ‘ 4 
B MN col)N nuajer megarayleighs/micron (2.4) 


where the resolution (AA) is given by 


Mi_ = 0.013 +(1.4 875 x 1077 A 


2 


0.0256 +(1.20 x lo “)A 


and 


~3 
= < 
HIGH GAIN Fy ) 6.5 x 10 Vw > Viry = 10 volts 
= 0.556 Vin 


» Vm, £1.8 volts 


F 
(Vig) me 


LOW GAIN 


= 0.sig9y°71277 


| F 


va < 
1.8 Vary © 20 volts 
The values of M 
oe ie 
Table 2.1-1. Mcol)N 
limits. The AFGL calibration was truncated based cn data at 2.10, 3.50 


(col)N depends upon wavelength and are given in 


are set to zero for values of A’ outside the tabled 


and 3.25, 5.40 microns. ‘These data are plotted in Figure 2.1-1 and show 

the wavelength overlap around 3.25 microns. This fact provided a conveni- 
ent region to check consistency in calibration between the long and short 
wavelength filter halves. Notice the brightness scale on Figure 2.1-1 is 


10° MR/um. Thus the calibration peak at 2.95 um is 8228 MR/ym. 


Telemetry Signal 


The EXCEDE II employed four telemetry (TM) links with ground 
receiving and recording systems. All data for the CVF instruments was 
contained in TM Link 3 tape (recorder 2 original) track S of the 
range tape. IRIG B time code to be used with these data was on track 
3 station MUX (the normal 13 IRIG B track was reported to be spotty.) 

The data was pulse code modulated (PCM) which is described by Space 

Data Corporation (1979). Data for each instrument coded in several 

PCM words was "read", reformatted and recorded on separate digital tapes 
along with time in a format compatible with the AFGL CDC 6600 computer. 

Flight data was digitized from T+40 sec (0457:20 UT hours) to 
T+340 sec (0553:50 UT) which was loss of signal. Data for the \S-6-4 CVF 
instrument were in words 2 and 26, 3 and 27, and 4 and 28 of all 
frames, Words 2 and 26 are high gain, 3 and 27 are low gain and 4 and 28 
are the reference. The high and low gain contain the spectra and the 
reference contains the start pulse of each spectral scan, in addition 


to other wavelength reference pulses (of lower amplitude than the 
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TABLE 2.1-1 
M(col),, vs. Wavelength NS 6-4 


Short Wavelength Half 


.03 x 102 
.03 
195 
C77 
.56 
.40 
.28 
.24 
27 
.38 
51 
.60 
.94 
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.01 
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43 
£15 
.81 
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.16 
91 
74 
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Long Wavelength Half 
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91 
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start pulse) that were not used in the data processing. A review of 
these words showed a one data point lag in the sturt pulses between 
words 4 and 27. To avoid this confusion only word 4 was used to 
locate the scan start times. 

Criteria for scan start was the first data point of value 20 
counts or less following the peak (?900 counts) of the start pulse. 
A listing of these start times (with other parameters is given in 
Appendix A). 

The high and low gain spectral data were contained in words 2 +26 
and 3+27 were converted to TM volts using the equation 


C-12 
pallets 2-8 
™ = 100 (2-8) 


Where C = number of counts (Space Data Corp. 1979). 


Data Processing 


The CDC digital tapes described in section 2.2 were processed 
using the scan start times (Appendix A) and the calibration equations 
to form a data base of one record per spectral scan. The seven 


parameters per data sample are 


Parameter Content 
1 wavelength - microns 
Z time after launch - seconds 
3 altitude - km 
4 high gain TM voltage - volts 
5 high gain brightness - MR/ym 
6 low gain TM voltage - volts 
7 low gain brightness - MR/ym 


These parameters were computed for every data point. Payload altitude 
was computed with an eighth order fit to radar range data and utilized 


time after launch as the independent variable. The trajectory equation 


and coefficients are 


° 


ee 


where 


i=8 Fl 
h=a,+ ££ a.T 
0 1 
i=l 
ay = -.539640506E+02 ae = - 1480095146E-07 
ay = .2741505348E+0)] ay = -.4493096965E-10 
a, = -.294716920SE-v1 a. = .7464427803E-15 
/ 
a, = .3475921337E-05 ay = -.5216526096E-16 
a, = -.2907943273E-05 


h = payload altitude - km 


T = time after launch - seconds 


The brightness data was plotted on semilog scale for each gain 
as a function of wavelength and altitude. All scans were plotted. 
The data was also plotted on a linear scale but otherwise in the same 
format. A sample brightness spectra is shown ir Figure 2.3-1. The 
data shows Wu Benesch spectra between 2 and 3 microns and the prominent 
4.28 micron cO,. In the 3.4 micron region the filter short and long 
half responses overlap and indicates good agrcement in the absolute 
brightness calibration. Notice that every spectral data point occurs 
at a different altitude because of payload motion as reported previously. 
Below 250 MR/um the presence of noise equivalent to 1 count in the 
onboard digitization process is evident. This effect is described 


in the next section of this report. 


2.3-1 Error Analysis 


For the NS-6 CVF the uncertainty in establishing counts is ] 


count, therefore from Equation (2-8) 
min AV = 0.01 volts 


Thus for high and low gain F(V,,,) is 6.5 x 10-5 and 5.56 x 1073 
respectively see Equations (2-6) and (2-7). For the short half the worst 
case uncertainty occurs at 2.95 micros (see Figure 2.1-1) and the 
uncertainty for high gain is 0.53 MR/um. For the long half the worst 

case high and low gain uncertainty occur at 3.25 microns and are .03 MR/um 
and 24,7 MR/um respectively. These values represent error bars {+ and 


-) of the measured data and in gencral are given by 


a F(.01) (2-10) 


where 


Fo) = factor for Vim = 0.01 volts 


Mecol) 


wl Figure 2,1-1. 


The percentage error represented by the one error count depends only on the 


number of counts measured, viz, 


AB 1 
- oP Ses 2, 2. 
€ x 100 fiz * 100 (2-11) 


It is desirable to merge the low gain data with the high gain data when 
the high gain saturates. However as will be shown below problems of accuracy 
are incurred. High and low gain volts are related throu .. the electrons gain 


factor G as follows 


V,, = GV (22725 


OQ 
" 


gain factor 


< 
1 


Ho high gain volts 


< 
" 


low gain volts 


From Iiquation (1) the nigh gain counts (C5) and low gain counts (Cy) are substi- 


tuted in Equation (10). 


Git 12061) 9G, 
Cae ee ag + de (2<15) 


If the switch from high gain to low gain is done using equation (11) the error 
due to 1 comt at the switch over to low gain is found by substituting (11) in 
(9) which yields 


: a x 100 (2-14) 
H 


“switch 
If we choose to switch at CG, = 900 the crror after the switch due to 1 count 
uncertainty is 8.3% assuming G = 75 where just prior to the switch the error is 
only 0.11% due to 1 count uncertainty. This discrepancy at the switchover point 
is undesirabie and consequently the high and low gain data was not merged but 


treated separately. 
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2.3.2 Data Bases 


In view of the incompatibility of high gain low gain data merge described 
above, no attempt wus made to merge the two channels. Instead high and low 
gain channels were used to generate individual data bases. Samples from these 
data bases are given in Appendix D. High gain channels show the effects of 
one count digitization noise when the guns (accelerators) are not in operation, 
superimposea on a background signal caused by a small offset voltage in the 
system (see scan 116, 119, 121, etc.). The spectral structure of the background 
mirrors the system responsivity characteristic (see Figure 2.1-1). 
In the region of 100 km (scans 139, 153, 159) there appears to be a substantial 
signal enhancement probably due to instrument covers floating through the CVF 


field of view. 


Periods for gun-on operation are given in Table 2.3-1. The correlation 


of signal enhancement with periods of gun-on operation can be observed by 


inspecting scans 237, 276, etc. 


Table 2.3-1 


ACCELERATOR OPERAT10ON 


Gun On Gun Off 
Time* Time* Gun 
(seconds) (seconds ) Number 
115.495 116.052 1 
135. 845 136.104 4 
138.651 139.695 4 
141,636 144.148 4 
153.254 156. 745 4 
159.277 161.547 4 
173.220 174.538 4 
176.485 176.728 4 
179.274 100.335 4 
; 188.113 191.242 4 
193.916 197.811 4 
208.083 209. 383 4 
211.355 212.549 4 
213.885 214.625 3 
222.984 226. 814 4 
228.779 232.681 4 
240.424 240.325 3 
243. 302 244.273 3 
243.032 244.285 4 
248.724 248.928 3 
240.246 250.159 4 
257.881 281.735 3 
ae 257.885 261.743 4 
283.607 263.925 1 
263. 704 264.331 4 
263.704 260.7 3 


Table 2.3-1 (Cont.) 


Gun On Gun Off 
Time* Time* Gun 
(seconds) (seconds) Number 
275.367 279.228 4 
281.083 282.047 1 
281.189 282.750 4 
: 234.525 235.032 1 
: 283.684 285.088 3 
286.929 290.859 l 
292.863 293.236 4 
292.759 244.156 1 
298.149 296.72 2 
292.863 296.723 3 
298.581 298.795 1 
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PART B 


He CVF SPECTROMETER 


3.0 Description of Helium CVF (HS-3B-1) 

The HS-3B-1 CVF (see specifications Table 3.0-1) is basically similar to 
the nitrogen CVF except a tuning fork optical chopper was added that periodi- 
cally interupted the incoming radiation to the sensor at a frequency of ap- 
proximately 106 Hz. The sensor electrical output was filtered with an elec- 
tronic filter cerrered at 100 Hz and passband approximately 50 Hz wide. This 
technique produce: an amplitude modulated waveform whose envelope contained 
the observed spectral radiance information. Both a high and low gain channel 
were included f¢r the modulated signal. The basic problem in processing the 
Helium CVF data was to demodulate and rectify the output signal then apply 
calibration factors to deduce the observed spectra in MR/um. 

Also unlike the Nitrogen CVF the long half of the HS-3B-1 CVF consisted 
of two filter segments in order to achieve the long wavelength response. How- 
ever since the long half wavelength calibration applied throughout the spectral 
region the brightness calibrations at the joint wavelengths was simply set to 


zero and the long half of the CVF was treated as if it were continuc. 


3.1 Calibration Parameters 


Radiance and wavelength calibrations were accomplished by both AFGL 
(Condron and Austin, private communications) and USU (Wyatt, 1979). Agreement 
of these data was adequate however the AFGL calibration data were used because 
the information was more detailed. In addition the AFGL wavelength data was 
also slightly modified (shifted) to agree with certain known emission features 


in the data such as co. 14.97 wm. However the shift was less than a wavelength 


Table 3.0-1 


PERFORMANCE PARAMETERS FOR CVF HS-3B-1 


Wavelength Coverage (um) 3.9 to 6.8 , 12.7 to 22.3 
Spectral Resolution (um) 0.16 @5.3 pm , 0.4 © 15 um 
Detector Type Si : (As) 
J s 2 * 

NESR (w-em72-sr Lim 1) 1.08 x 10 i] 

: 4* 
Dynamic Range 8 x 10 
Half Angle Field of View (deg) 1.973 
Spectral Scan Rate ~0.8/sec 


* 
reference 7 


resolution element of the instrument, 


3.1.1 Wavelength Calibration 
The HS-3B-1 filter halves are designated long wavelength (subscript ¢) 
and short wavelength (subscript 4). The short half was scanned first (in 


time}. The adjusted calibration equations for wavelength for each half are 


given by 
dp = (2.388 E-1) Pa - 1.097 um (3-1) 
55 ¢ P, s 97% 
d, = (7.38 E-2) = + 3.333 um (3-2) 


where Es is the percentage into the scan ae = 100 n/N). 
The wavelength resolution elements at any wavelength determined by 


the CVF filter are given for each filter half by, 
AA, = (1.8667 E~2) A + 6.21E-2 um (5-3) 


Ap = (5.482 E-3} } + 3.12E-1 um 


3.1.2 Radiance Calibration 
During actual radiance (brightness) calibration of the HS-3B-1 CVF the 
: sensor output was a continuous analog signal and demodulation and rectifi- 
cation of the signal was accomplished using a phase sensitive amplifier 
(PAR). Since the telemetry signal was PCM this technique was not possible 
so a computerized demodulation and rectification scheme had to be formula- 
ted to accomplish these tasks which included the laboratory calibration data. 


This formulation is derived in section 3.3.1 of this report. However, the 
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final result to convert demodulated rectified telemetry volts Vom to spec- 


tral radiance is given here 


By = G; M(co1)“ry/o4 MR/pm (3-4) 


where BL brightness (radiance) of data point n 
G; = channel gain (high gain = 1, low gair = 100) 
M601) = calibration factor in MR/volt (see Table 2.1-1) 


Ad = resolution (see equation (3-3)) 


The adjusted laboratory calibration data for M (col) is given in Table 3.1-1 
in units of MR/volt. Notice that the values are zero at 17.50 and 17.70 um 
where the long wavelength filter point occurs. Other values in the original 
data near the mask regions are located have also been adjusted to eliminate 
the tendency to "spike" when the mask passes over the sensor. A plot of 

the final calibration is shown in Figure 3.1-1 in units of MR/um - volt and 


derived using the data in Table 3.1-1 and equation G-3. 


3.2 Telemetry Signal 


Beg All aspects of the HS-3B-1 PCM telemetry were identical with CVF 
NS-6-4 (see section 2.2 of this report) except that since the HS-3B-1 CVF 
was optically chopped one additional channel of data was required (compared 

: to the NS CVF) for the chopping frequency. The spectral data was contained 
in PCM words 5 + 29 (low gain) and words 6 + 30 (high gain). The spectral 
reference (start pulse) was on words 7 + 31 and the chopper sequel on word 
46. Spectral scan start times are determined in the same manner as the NS 


instrument, see Appendix B for the listing of HS-3B-1 start times. 
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3.70 
3.80 
3.90 
4.00 
4.10 
4.20 
4.30 
4.40 
4.50 
4.60 
4.70 
4.80 
4.90 
5.00 
5.10 
5.20 
5.30 
5.40 
5.50 
5.60 
5.70 
5.80 
5.90 
: 6.00 
6.10 
6.20 
6.30 
6.40 
6.50 
6.60 
6.70 
6.80 
6.90 


1 
M(col), 
MR/volt 


Short Wavelength Half 
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TABLE 


3.1-1 


M(col),, vs. Wavelength HS 3bB-1 
M(col),, » 
MR/volt pm 


X 


um 


Long Wavelength Half 


0 18.50 
0 18.70 
0 18.90 
0) 19.10 
53 19.30 
26 19.50 
.16 19.70 
15 19,90 
.09 20.10 
.0S 20.30 
.025 20.50 
.9635 20.70 
. 8545 20.90 
.§36 21.10 
-8855 21.30 
. 8905 21.50 
.835 21.70 
.7925 21.90 
. 7665 22.10 
.763 22.30 
. 804 22.50 
.853 22.70 
8765S 22.90 
8205 
716 
.735 
7545 
8725 


M(col),. 
MR/volt 


8375 
8225 
82) 
843 
. 869 
8825 
902 
9235 
8985 
875 
8935 
920 
9405 
£951 
964 
9805 
9805 
9805 
£9805 
9805 

0 

n) 

0 
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Data Processing 


The idea of chopped systems accompanied by phase sensitive demodula- 
tion and detection to reduce DC drift and improve signal to noise is well 
known (Schwartz, 1959) and widely used by electronic engineers. However 
for actual digital CVF data such as that of NS-3B-1 sore special problems 
occur and so a brief review of the fundamentals wil! be presented prior 
to describing the actual data reduction techniques used for CVF HS-3B-1. 

Figure 3.3-1 illustrates how the HS-3B-1 CVF modulated signal was pro- 
duced. Since the CVF selectively passes the incident spectral radiance as 
it rotates at a constant rate the incident spectra! radiance can be repre- 
sented as a function of time instead of wavelength. This is shown in the 
top plate of Figure 3.3-1. The incident spectral radiance is periodically 
interupted by a mechanical chopper which itself radiates some lower level 
since it is cooled. The chopped signal as seen by the CVF sensor is shown 
in the second plate of Figure 3.3-1. The sensor output is assumed propor- 
tional to this incident radiance. If we assume for the moment that the 
chopper radiance is negligible and that the frequency bandwidth of the in- 
cident spectral radiance is small compared to the chopper frequency then 


the sensor output can be estimated by the Fourier series 


s(t) =F S,) [ + 2 2 Sin le cos nut] (3-5) 


chopped spectra signal 


where s(t) 


yn 
il 


(X) infrared spectra (unchopped) 


radian chopping frequency 


€ 
i] 


If we now assume that there is no incident spectral signal and the 


chopper radiance is not negligible then the signal output from the chopper 
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radiance alone can be written in a similar manner except that there will 


be a 180° phase difference 


ao 


_1 Sin n 7/2 - 
es) = 7 Ry) ft + 2 » ies. COS (nagt + 7) | 


n=1 
where Dee) = chopped chopper signal 
Roa) = chopper unchopped radiance 
uF radian chopper frequency 


Clearly the phase angle 1m can be omitted simply by changing the + sign in 
front of the summation to -. 
For the case of both incident spectral radiance and non zero chopper 


radiance the total sensor output signal is the sum of s(t) and bee) 
f~ co 


g(t}=8(t)#b(_)25(Sqy*Raay)*Sy-Ryay) > ane cos nut (3-6) 


n=) 
If we now A.C. Couple this output signal to remove the DC term #54) + Riay? 
then filter with a narrow bandpass filter centered 4. the fundamental chop- 


ping frequency wy we get for n=] 
(t) = = (S - R,,,) cos wt (3-7) 
: 7A)” *Q) c : 


This result is the modulated waveform of the EXCEDE sensor output shown in 

plate 3 of Figure 3.3-1. This sensor output was digitized on the payload and 

telemetered and recorded in PCM format. Actual data is shown in Figure 3.35-1A. 
The sinusoidal waveform of the chopper reference (also digitized and 


telenetered) in plate 4 of Figure 3.3-1] is simply expressed as 


h(t) = A cos (wot + Aé@) 


where h(t) = A.C. Component of the chopper signal. 


A® = phase difference between data and chopper signals. 


A = amplitude of the chopper signal. 
The phase difference A@ is included in this equation to account for such 
things as mechanical lag in sensing chopper position which would produce 


a phase difference between the sensor signal and the chepper reference. 


In principle the demodulation and extraction of the spectra S 


(A) 
from equation G-7) is done by forming the product of h(t) and g(t), which 
expanded and simplified yields 


7 (ee ee 
P = g(t)h(t) = ———O)_ tr)" cos wt cos (wt + A@) 


(3-9) 


A (S,.,-R,,,) 
= — [cosas + cos(2u t + ae)] . 


We now filter the product P with a low pass filter with cutoff frequency 
less than 20, and greater than the highest data frequency in the spectral 
radiance. Solving for S(, yields 


1P* 


= Roy) * Kos be 


(3-10) 


where P* represents the low pass filtered product of h(t) and g(t). 

Equation 6-10) shows that the chopper radiance background is added to 
the incident signal and clarifies why it is necessary to keep the chopper 
radiance small (cold). Also if either the chopper amplitude A or the 
phase angle difference A€@ between spectral data and chopper reference are 
modulated with time in any manner false components will be introduced into 
the reduced spectra. 

Notice from Equation 6-7 that if the incoming radiation Sian) is less 


than the chopper radiation Rey) then g(t) is negative and the product term 


"a. a ne lle $e aaa 


in Equation G-10) will also be negative resulting in a quantity of radiation 
being subtracted from Reay to give re) just as it should be. If one 


attempted to retrieve the spectra S by simply reconstructing the envelope 


(A) 
of g(t) this possibility would not be detected and an error in 314) would 
result. On the other hand multiplying by chopper signal and filtering as 
explained, preserves all the phase properties of the signal. It is crucial 


however, to derive an accurate representation of the chopper data character- 


istics. 


3.3.1 Data Reduction 


Several problems arise in attempting to apply the principles derived 
above to the EXCEDE chopped CVF PCM data. For one thing all the data 
channels were sampled (on the payload) serially so that the chopper refer- 
ence data points were not time coincident with the spectral data points. 

In addition the chopper was sampled at 512 Hz - just half the spectral data 
sample rate. Thus forming the correct product P (Equation (3-9}) of the 

chopper and spectral data needed for demodulation was not possible. Con- 
sequently it was elected to fit the digital chopper data for each spectral 
scan with a least squares technique to derive an analog equation to be used 

in rorming the product term in Equation G-3. The mean value of the chopper 
reterence signal data was found and subtracted from the data set. An estimated 
frequency (f.) of the signal was found by counting cero crossings during 

each scan and dividing by two times the product of the total number of data 


points and the sampling interval, viz. 


NAf 


Q 0's 
fo = > = (3-11) 
c Neg ht 2Noor 


ty 
tu 
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where fo = chopper frequency - Hz 
No = number of zero crossing in the scan 
Nror = total number of data points in the scan (1/2 number Spectral data 


samples) 


At = sampling interval 


f_ = 1/At = sampling frequency = 512 Hz. 


Since the total number of chopper data points Neor 


the number of spectral data points) for each scan and the sampling frequency 


is known precisely (1/2 


is considered invariant the error Af, in determining f. is at most due to 


the error of missing one zero crossing and thus 


_ (1)512 


Af. see) * 0.42 Hz 


(The value 616 for N, was estimated based on 1232 data points in a typical 


TOT 
spectral data scan - see Appendix A). 
Knowing the frequency the assumed chopper sinusoidal signal equation 


derived below was least squares fitted to the data in each spectral scan. 


h(t) 


Acos (wt -9. )= Acoswt cos8 + A sinw.t sind. (3-12) 


h(t) = b cos wt + ¢ sin wet 


where h(t) = chopper signal 


A = amplitude 
; ead phase angle relative to the scan start 
b = 
b A cos o. 
¢ = A sin 6, 


As a result of the fit the constants b and c are found and the phase 


angle is determined by (see Figure 3.3-2 for typical data fit) 


8 = stan) cS) (3-13) 


a’ 
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Finally setting the amplitude A=1 the chopper signal equation in digital 


form is 


where h(n) = chopper equation 


ft. = chopper frequency - Hz 

n = data point number 

At = sample interval = rosy sec: for the spectral data 
8. = phase angle relative to the scan start 

Noor = total data points in the scan. 


In actual practice tne frequency a was computed by the fit. Appendix C 
contains a tabulation of the pnase angle and f. ( for the actual data derived 
as discussed above for each spectral scan. Included in this tabulation is 
the "goodness of fit" parameter chi squared. The lower the number the better 
the fit. <A typical example of the calculated chopper reference (solid line) 
using Equation 6-14) compared to the actual chopper datu (points) is shown in 
Figure 35.3-2. As iong as the chopper reference data umplitude was relatively 
constant (which was not always the case) the fit with Equation 6-14 was 
excellent. 

Clearly Equation G-i4 permits calculating the necessary value of the 
chopper reterence for any spectrai data point n and thus satisfies the re> 
quirements needed for evaluating the product P in Equation @-9) (which impiied 
579). When there is a phase difference between the spectral data and the 
chopper defined by ire = £0 where Ms and U, are the phase differences from 


the start of scan for spectra and chopper respectively, Equation (5-7) becomes 


g(t) = FO oy ay cos(w t-C.) 


It was found by inspection of the data that the chopper lagged the spectral 


data by 1/3 a spectral data point (12°) so that A9 = -12°. 


An example of the product P for scan 1211 formed by multiplying each 
spectral data point in the scan by a numerical value computed for the chop- 
per using Equation 3-14 for the same n (constant A=] and the appropriate 
phase angle e., taken from Appendix C) is shown in Figure 3.3-3, both the high 
and low gain channels are shown. The product was evaluated in telemer-ry 
units of counts for the spectral data g(t) while h(t) is dimensionless. 
Prior to forming the product the mean of the modulated spectral data was 
determined and subtracted from each data value of g(t). Theoretically the 
mean of the modulated wave should be zero. Note that all prominent count 
values, in Figure 3.3-3, are positive except for a few spurious data points 
in the high gain channel where the signal was saturated. In addition the 
noise values oscillate about zero as it should. The envelope of the data 
in Figure 3,.3-3 contains the desired spectral information. All that is re- 


quired to derive the spectrum is to extract the envelope and apply calibra- 


tions for wavelength and brightness. Ripple on the high gain is explained later. 


Typically envelope detection is done by smoothing or filtering. 
Clearly the only useful data points of the waveforms shown in Figure 3.5-3 
are the peak values since only they represent the observed spectrum. Since 
filtering of digital data is a somewhat complicated process especially with 
large data bases and when Fourier transforms are involved it was decided 
to attempt to derive the envelope data by a decimation technique prior to 
filtering. | 

Using the concept of referring phase angles of both the chopper reference 
and spectral data to the beginning of the scan and acknowledging that these 
phase angles could be different the product term of Equation G-9) can he 


rewritten a5 


| P= AS Gy -Reay) eos (85-8) + cos [2u,t-(8 +.) 1f (3-15) 


since eae = A6 is constant the maximum value of the product P occurs when 
the cos [2w.t-(6.+8.)] = +1. Thus the data point number ny defining the maximum 
value is 
2{2z(n,-1) fF At] - (6.+95) = 2(1-1) 7 
121525350415) GNaAas oo 


Therefore solving for a, yields 
(i-1) 7 + (6 .+6.)/2 
pa ash (3-16) 
p 2uf At : 
¢c 
Since Le is arbitrary it is clear that values of ny computed with Equation 6-1q 
- . could take on non-integer values which do not exist in the digital data base, 


viz. data points ure numbered in integer steps from 1 at scan start to N,... at 


TO? 
the end of the scan. Assuming the sampling rate fe is invariant (1024 samples/ 
sec) and equal to the reciprocal of the sample interval time At Equation 6-1) 


is rewritten as follows 


(3-17) 


cr (1-1) 7 + ee] 
o 


n= round oft integer value |1 + mise 
ies ~ 


a i = 1,2,3,...Nio¢ . 
The envelope data base was determined for each scan by choosing only data 
. points defined with Equation @-17. The data point difference An, = 5 and thus 
moos the data buse was reduced by a factor of 5 by this decimation process and con- 
tained only the estimated envelope data. 
Once the envelope data base was generated in this manner the spectral 


data in counts was converted to volts using the relationship 


ae ey 4 
Vom * Too * Too (e748) 
where Vow = telemetry volts 
C = counts 
26 
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This equation differs slightly from that used for the NS instrument (Fquation 
(2-8)) because the A.C. coupling of the HS sensor output removes the 12 count 
D.C. offset. The decimated envelope data base derived by Equation G-10 was 
units of counts. These data were converted to volts simply by dividing by 109 
in accordance with Equation (@-19. 

Returning to the basic relationship for deriving t':e product p in Equation 
(3-19, we must recall in subsequent discussion that we have already assumed 
the amplitude A of the chopper reference signa! to be constant and equal to } 
and also that the phase difference A between chopper refcrence and srectral 
data to be +12°. Finally we make the assumption that chopper radiance AoH 
is negligible compared to the spectrum since it is cooled. In addition we 


have used ny in choosing the envelope instead of Le so solving for the spectrum 


Tey) yields 
Sears ud (3-19) 
~ Sosh® + costdlonf (n -l)AtT- 1" -© V/J21 10 
(A) ~ cosh® + cos2[2nf (nj-1)At - (%. df 2] 
The units of 0) are counts. Notice the denominator -2 for small 4°. It was 


decided to first generate a data base where the spectra was not filtered but 
in the correct units of MR/um then filter the data set in units ef radiance. 
This precaution was taken in the event the voltage calibration was not linear. 


The unfiltered spectral radiance was thus calculated from the unfiltered 


Sa) 
decimated data base where Vom = Too 
} ’ a 
Ea Beay = 55 Mecor) View? (3-20) 
where Boy) = unfiltered spectral radiance MR/ym 
| P = decimated product term (unfiltered) in units of 
leg. = cultbration term - MR/volt 
A\ = calibration term - um (see Equation G-3)) 
Vim = decimated spectra term expressed in volts 
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G. = gain factor (high gain = 1, low gain = 100) 


Equation 6-21) was given previously in Equation G-4) of this report. 

Figure 3.3-4 shows a sample unfiltered spectral radiance profile derived 
by the techniques described above. The ripple on the spectra is caused by 
a chopping feed-through effect (viz cos 2u.t term in Equation 6-9)) which was 
intentionally not filtered before initial examination of the spectral data. 
The spectral data shows the emission features of co, at 4.28 and 14.97 ym 
and NO in the 5 to 6 um region. However the need of filtering to remove 
the chopper effect is evident. 

As a first approach a simple nonrecursive filter was used for initial 


data smoothing. The filtered brightness Be at point n is computed as 


follows 


Bos [aB > + bB 4 + cB + DB yy + aB 42! (3-21) 


The transfer function of this filter is (Hamming, 1979) 
How) = 2a cos 2u + 2b coswte (3-22) 
where w = 7 f/f, 
fy = 1/2 At = fs/2 Nyquist frequency 


fs = sampling frequency after decimating = 205 Hz 


f = physical frequency 


For this case we desire good low pass characteristics, so using the condi- 


tions 
H(O) = 1 
dH (0 ey ee 
es Sonn 
d7H(0) _ 
7) = 0 
dw 


leads to the following values 


a= -1/16 
b= +1/4 
c = +5/8 (3-24 ) 


A plot of the frequency response of this filter is shown in Figure 3.3-5. 

The decimated data base was processed using Equation 6-21) to produce 
initial smoothed data. Spectral scan 1210 (Figure 3.3-6) illustrates the 
effect of the smoothing and the improvement of spectral radiance features 


are shown in Figure 3.3-6 as compared to Figure 3.3-4. 


3.3.2 Data Base (HS~3B-1) 

An "abridged" data base was generated for selected spectral scans from 
the HS-3B-1 CVF using the procedures outlined above. These data are given 
in Appendix E. In general the spectra show emission enhancement during 
periods of "gun" (accelerator) on times but in addition there seems to be a 
background enhancement (scan 1067) at times when the guns are off. This is 
attributed to the possibility of instrument covers floating through the CVF 
field of view. The fine structure in some of the enhanced emissions (see 
5.8 um of scan 1204) is due in part to the effects of chopper second harmonic 
and sampling which is explained in more detail in the next section of this 


report. 


It should be noted that the spectral data shown in Appendix E consist 
of merged results from both the high and low gain channels. A signal voltage 


of 4.2 volts was used as the criterion for merging the two gain channels. 
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3.4 Discussion 


The ripple (~40 Hz) supperimposed on the demodulated waveform in 
Figure 3.3-3 is clearly an error source since it falls within the digital 
fiiter passband (see Figure 3.3-5). Power spectral density analyses of these 
Signals do not indicate the presence of the 40 Hz ripple, however, the 
chopper second harmonic (212 Hz) is evident in the PSDs of both the spectral 
and chopper data. This indicates that the chopper second harmonic is leaking 
into the signals due to poor frequency cutoff characteristics of the system 
electronic filter on-board the payload. Inspection of the electronic charac- 
teristic in Figure 3.4-1 confirms this possibility. A simple trigonometric 
exercise demonstrates the point. Assume tne chopper fundamental waveform 
is added to the second harmonic waveform approximately 90° out of phase and 


with reduced amplitude "a". 


h = COSw t + asin(2u t-s) (3-35) 


(t) 

where 9 = small phase angle -9° 

h = chopper wave form 
f= chopper frequency 106 Hz 
w = 2n(n-1) f/f, 
f. = sampling frequency 512 ilz 
n = data point number 
a = amplitude <1 


0.2 and evaluating Equation (3-39 for n=54 to 154 we produce the 


é Settiny a 
Simulated chopper waveform in Figure 3,4-2 (scale factors have been introduced 
for ease of comparison). Comparing this to a portion of the actual waveform 
indicates obvious similarity. Note that the ripple frequency in the chopper 


is only ~20 Hz (rather than 40 Hz as it is in the data) because sampling of the 
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chopper was only 512 Hz compared to 1024 Hz in the case of the spectral data. 

As pointed out the ripple cannot be removed from the data without more sophis- 

ticated filtering techniques which are beyond the scope of this preliminary work. 
There is also an inherent error in the spectral data due to a one 

count uncertainty due to digitization. One count is equivalent to 10 milli- 

volts by Equation 3-18. From Equation 3-20 the noise radiance for this 

voltage in high gain (G=1) is 


Bain = 9-91 Mecgyy / OA (3-36) 


where Mecot) / Ad is the calibration factor in MR/volt (see Figure 3.1-1) and 
Ak is the resolution. From the figure the calibration factor depends on 
wavelength showing values of 2.23 MR/um-volt at ~15 um and 16.9 MR/um-volt 


at ~5.3 um. This produces a noise of 


Ds 0.02 MR/um at 15 ym 


Bie 1.7 MR/um at 5.3 um 


These values are well below the signal levels encountered (see Appendix F). 
An additional error is introduced by the decimation process because of 

the necessity to round off the computed peak value data point number to the 

integer value. This causes a cyclical error every S data points of about 


10%. 
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PART C 


4.0 RECOMMENDATIONS 

Data processing techniques of Circular Variable Filter Spectrometers 
have developed in the past 5 years to a point where spectral resuits can be 
accepted with good confidence. However, limitations in interpreting 
processed spectra are imposed by a variety of design parameters. For instance 
leakage due to imperfect out-of-band rejection of the CVF can contaminate iow 
spectral signals. The best blocking possible with a CVF is probably four 
orders of magnitude. Non-constant CVI’ rotation can cause wavelength errors 
comparable with the spectral resolution of the CVF. In the case of the EXCEDE 
HS-3B-1 poor rejection properties of the electronic filter caused chopper 
second harmonic to contaminate the data. 

Specifically we recommend that the EXCEDE data and techniques presented 
here be considered preliminary and reviewed for adequacy in terms of the 
physics involved, The helium instrument (HS-3B-1) data should be reprocessed 
by more sophisticated techniques to reduce as much as possible the effects 


of second harmonic contamination. 
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Figure 2,1-1, CVF brightness calibration (NS-6-4) 
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Filter Response - H(w) 


Figure 3.3-5,. 
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Frequency response of smoothing filter used 
after decimating data base to a sampling 
frequency of 205 Hz EXCEDE information band- 
width is 50 Iz. 
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Figure 3.4-1 Characteristics of the electronic 
filter for CVF HS-SB-1 
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Figure 3.4-2. Comparison of chopper data with model of the form 
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Other Data for NS-6-4 CVF 


Data in this appendix is defined as follows: 


3 Column Title Definition 
i: 1 Scan Number Number assigned to the spectral scan 
I 
2 Length Number of data points in spectral scan 
a 3 Length sec Duration of the spectral scan in 
seconds 
. 4 Time A-.L Time after launch in seconds of scan 
: start 

5 Altitude Payload altitude corresponding to 

T.A.L. 
6 Time U.T. T.A.L. in U.T. seconds 
7 Clock Time U.T. T.A.L. in clock form 
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APPENDIX B 


Spectral Scan Start Times and 
Other Data for HS-3B-1 CVF 


Data in this appendix is defined as follows: 


: Column Title 

; 1 Scan Number 

fu Length 

rt : 3 Length sec 

4 Time A-L 
5 Altitude 
6 Time U.T. 
7 Clock Time U.T. 


Definition 
Number assigned to the spectral scan 
Number of data points in spectral scan 


Duration of the spectral scan in 
seconds 


Time after launch in seconds of scan 
start 


Payload altitude corresponding to 
TsAvLe 


T.A.L. in U.T. seconds 


T.A.L. in clock form 
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APPENDIX C 


Chopper Reference Phase Angle and 
Other Data for HUS-3B-1 CVF 


Constants A, b, and c are defined by Equation 3-12 of the text. The 
frequency (FREQ) is the best fit chopper frequency and THETA is the chopper 
phase angle O., relative to the scan start. MEAN HI and LO are the spectral 
mean counts for the high and low gain. NPTS is the number of points of the 
scan. CHISQ and INTER are the yoodness of fit parameter and the number of 


interations required to fit the data. 
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APPENDIX D 


Sample Data Base for CVF NS-6-4 


This appendix contains selected spectral data developed from CVF NS-6-4 


measurements. Scans and scan start times (time after launch in seconds) It 


follow 
Scan TLALL. Scan Teele 

116 96. 1080 271 ~—-:171.4120 
: 119 97.5590 276 «173. 8840 
: 121 98.5300 395 232.7120 
ee 126 100.9500 424 247.0380 
136-105. 7810 447.257.6770 
139: 107.2340 467 207.3540 ' 
iS3 114.0030 487 277.0310 
156 115.4520 -497 281. 8680 
159 116.9030 502 284.2890 
171 122.7050 505 285.7450 
174. —-:124.1570 507-286. 7110 
181 127.5440 510 288.1610 
223. «147. 8600 -512 289.1290 
228 150.2780 520. 292.9960 } 
231 151. 7300 2g02 293.965 


- indicates difficulty in detecting scan start time. 
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Figure 3.3-2. Example of least squares fit (sinusoid) to chopper signa: 
(solid squares; for spectral scan 1210. Also shown are 
frequency and constants determined as a result of the best 
fit. 
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Figure 3.3-5. Frequency response of smoothing filter used 
after decimating data base to a sampling 
frequency of 205 Hz EXCEDE information dand- 
width is 50 Hz. 
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